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ABSTRACT 


The critical stress intensity factors and fracture 


trajectories for both filled and unfilled viscoelastic plane 


stress specimens were experimentally deternined by 
controlled crack elongation. Fracture testing included 
combined displacement loadings both of Single and 


bi-material specimens with initial cracks located in each 


material and at the bi-material interface. 


The feasibility of trajectory prediction using an 
elastic linear-Strain finite-element analysis was 
investigated. Techniques for adapting the elastic model 
grid patterns to predict viscoelastic trajectories were 


developed. 


Measured trajectories correlated well with the 
directions of maximum principal stress as determined by the 


finite-element solution of a blunt crack-tip model. 
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I. INTRODUCTION 


A. BACKGROUND 


During the manufacture of solid propellant rocket 
motors, it is possible for stresses to develop during cure 
cooldown of sufficient magnitude to induce fracture at 
physical discontinuities and internal flaws. Such cracks 
may be tolerable providing the trajectories do not create 
Significant additional burning surface nor expose the case 
to combustion temperatures by local failure of the 
insulation. The direction of crack propagation is therefore 
of considerable significance. The stress levels for 
fracture propagation are of equal interest. It is therefore 
desirable to study the behavior of these viscoelastic 
materials and to develop technigues for evaluating their 


critical fracture criteria. 


The fracture stresses are commonly expressed in terms 


of the stress intensity factors K , K and K where the 
I Il Iil 


subscript refers to the fracture mode as shown in Figure 1. 
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B. CLASSICAL SOLUTION 


The study of fracture mechanics, to a large extent, has 
dealt primarily with the opening mode behavior of elastic 


materials. 


The classical solution to the detailed stress field of 
a sharp-tipped crack as shown in Figure 2 was formulated by 


Williams and given by Sih and Liebowitz [1] as: 


GD 


Figure 2. Stress Field Coordinate Systen. 


Symmetric Case (Mode I) 
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Y2r 2 2 Z (1b) 
oxy = KI cose[siné:cos36] +... 
Y2r 2 2 2 (1c) 
and 
u = ‘Var Ni Zed ieosil = Gosc6 | ot eee (2a) 
8G 2 Z 
v = “iv@r [(2x+1)sine - singe] +... (2b) 
8G 2 iz 


Skew-Symmetric Case (Mode IT) 





oe a Sia Sim @ [2G s Gre Os cule ae ee (3a) 
v2r 2 2 Z 
Oy = ATI Sabie eel) OOS Sell ae ao an (3b) 
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Oxy = Ki Eos | & jes mies nso + a.c (3c) 
Vor 2 2 2 
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and 





u = SiivZ [(2k+3)sin® + singe] +... (4a) 
8G is 2 
v = ATIV2r [(2«-3)cos@ + cos30] +... (4b) 
-8G Z Z 
where K = 3 - Av Paaners &raan 
= s oe = Plane Stress 
aes, 


George Irwin related the stress intensity factor to the 
strain energy release rate by computing the work done at a 


propagating crack tip. His results for K_ show that: 
L 


K.7= =f edU (5) 


By means of the energy balance introduced by A. A. 
Griffith in the early 1920's, it can be shown that res the 


value of K at fracture, must be a material property. 
I 


ip 





C. COMBINED LOADING 


Previous studies [2] concerning stress fields in rocket 
motor castings indicate eee equal shear and 
tension at the propellant-liner interface near the motor end 
Surfaces. It is therefore desirable to investigate the 
behavior of such materials under a combined loading 


condition. 


Adaptation of the classical solutions for single-mode 
stress fields to the combined-mode case requires only 
Superposition of the respective solutions. Combining the 


opening and shearing modes, we then have: 


z.. Sel cos@{K,[1-sin@+sin30]-K,,tan@[2+cos6-cos3e]} (6a) 
u2Yr 2 Z Z Z Z 2 

as ipcosoik, [I+sin8-sin36]*ky;tane[eose-coss0]] (6b) 
v2r 2 Z Z 2 2 Z 

oxy = 1 eae od (Kp i(sameicos sel Bass mick ce scre (6c) 
vcr Z 2 2 2 Z 


In polar coordinates, the stresses are given as: 


o. = 1 cos8[(3-cos6)K, a (3cos@-1)K,,tano] Hee a oa (7a) 
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Og = 1 cos6[(1+cos6)K, See) Sern x ay ae (7b) 








Oy ae 2 
Trg = BeOS) S10 mens cOs.0 hea eee (7c) 
DE 2 


It can be seen that the general formulation of the 


combined loading solution may be expressed as: 


Two commonly recognized hypotheses for crack extension 
in a brittle (non-viscoelastic) material under Slowly 
applied loads are: 

a) The crack extension starts at its tip in the 


radial direction, and : 


b) The crack extension starts in the plane 
perpendicular to the direction of greatest 
tension. 


These state that a crack will extend radially in the 
direction of maximum tangential stress where the shear 
stress is zero. Consistent with the theory of Griffith, the 
direction of crack propagation corresponds to the direction 


of maximum energy release. This direction has been 
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identified as the maximum principal stress direction. 


It 1S intended herein to demonstrate that the fracture 
trajectories of viscoelastic materials may be predicted by 
the maximum principal stress directions as determined by 
elastic analyses. TO Minimize the time-dependent 
viscoelastic effects, the loads were slowly applied in order 


to approximate an "equilibrium" condition. 


In order to determine the direction of maximun 
principal stress, the detailed stress field in the vicinity 
of the crack tip must be known. Williams' solution, based 
On an infinite sheet with a sharp-tipped central crack, has 
generally been accepted for any geometry in the "vicinity" 
of the crack tip. Definition of this "vicinity" varies with 
Sample geometry. Another method of establishing the stress 


- 


field is by finite-element techniques. 


There are two basic objectives of this study: 1) to 
demonstrate that the linear-Strain finite-element solution 
1s suitable for fracture analyses by correlation with the 
classical field for a sharp-tipped crack, and 2) to apply 
the finite-element technigue to determine the direction of 
maximum principal stress for the blunt-tipped crack 
geometry, which models the viscoelastic materials at 


fmeacture. 
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The results of crack elongation tests under combined 
loading with constant boundary displacements have been 
correlated herein with the trajectory predictions made by 


finite-element methods. 
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If. EXPERIMENTAL STUDIES 


A series of combined loading tests was performed on 
selected taterial specimens. The geometry of the specimens, 
whose fracture behavior is independent of crack length, is 
Shown in Figure 3. An initial crack length of 1.0 inch was 


chosen for all tests. 


A. DESCRIPTION OF SAMPLES 


Three basic specimens of the same geometry but 
different material properties were examined in this selane 
1) an unfilied viscoelastic liner material, 2) a filled 
viscoelastic propellant sample, and 3) a composite specimen 
which models the propellant-liner interface. Initial cracks 
were made both along the centerline and offset 0.2-inches 
from the centerline. The specimens were subjected to 
combined loading by constant boundary displacements ranging 


from pure opening (Mode I) to pure shear (Mode ITI). 


B. SAMPLE FABRICATION 


The Samples were fabricated by United Technology 
Center, Sunnyvale, California, The propellant and/or liner 


materials were cast in redwood boxes. The wood, which was 
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Geometry of Test Specimens. 





easily cut and shaped to the desired geometry, provided a 
semi-rigid boundary to the rubbery materials £OL 
installation in the testing apparatus. The inner surfaces 
of the wood were coated with an epoxy cement which 
penetrated the wood and provided a bonding agent for the 
liner. Since the propellant does not adhere to the epexy or 
to the wood, a thin film of liner (to which the propellant 
bonds well) was applied to the epoxy permeated surfaces. 
The samples were cured in the boxes and later sawed and 


milled into 0.10-inch thick specimens. 


The propellant samples were determined to be relatively 
uniform but the liner samples SOR SCeEea a "soft" region 
extending from the wood boundaries approximately 0.15-inches 
into the sample itself, resulting from incomplete cure near 
the extremities. The region was identified by a tacky 
surface with a glossy appearance. This undesirable feature 
‘resulted in severe necking, which precluded accurate 
determination of strain levels during testing. The 
effective gage length was reduced from the specimen height 
of 1.0 inch to a maximum of 0.7 inch for the single-material 
specimens and to a maximum of 0.35 inch for the composite 


Samples (on the liner half). The resultant force along the 
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direction of displacement was measured and modulus data were 


obtained from separate tests of the biaxial specimens. 


C. SAMPLE PREPARATION FOR FRACTURE TESTING 


The initial cracks were made by cutting a one-inch slit 
in the samples with a knife. The crack tip location was 
marked by grid lines drawn on the Sample itself with a fine 
point drawing pen (No. 00O Rapidograph). The grid lines 
facilitated detection of crack propagation and were used to 
determine the strain in modulus evaluation tests. The 
displacements of the grid nodes were optically measured with 
a cathetometer, since the "soft" areas invalidated 
cross-head travel data as a method of strain determination. 
Mechanical testing indicated the moduli varied by as much as 
50-percent between batches, and an equal variation within 
the same batch could be expected iQue different 
configurations and/or material aging time between tests. 
(The propellant is extremely hygroscopic. During initial 
tests the samples were retained in a dry nitrogen 
atmosphere, but due to the unavailability of a controlled 
environment during tests at NPS, the samples were allowed to 
age unrestrained. This affected the critical load eee 


but the trajectories are independent of absolute moduli). 
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D. LABORATORY PROCEDURES 


a Description of Hardware 


OF oer ={a—, er ee RF == ae ee a 


The combined-loading tests required construction 
of a specimen-holding apparatus, shown in Figure 4, which 
would allow rigid-body displacements of the specimen 


boundaries ranging from y = 0° to ye ee 





AANA AN ARNT 


Figure 4. Drawing of Specimen Loading Fixture. 
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2. Force Measurement Under Combined Loading 


ee ee a ee a omy eee ee ee ee ae ae eee eee ee 


Using the finite-element solution for stress 
distribution along the So aia akaa pavennal boundaries, the 
loading condition of a Sample was determined by measuring 
the force component in the direction of displacement. This 
force was transmitted to the testing machine load cell, and 
recorded, to establish the stress intensity factors. BY 
measuring the applied load vs. displacements, the influence 
of the "soft" regions at the liner-to-wood interface was 


Minimized. The loading diagram is shown in Figure 5. 





Pa . 
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Figure 5. Loading Diagram of Material Sample. 





3. Detection of Crack Elongation 


we a ae ee =o aw a Se —_ 2 <a er ee 


Identification of fracture loads waS hindered by 
the extreme low crack propagation velocities associated with 
viscoelastic fracture at low loading rates. During all 
fracture ‘tests the loading rate was maintained at 0.05 inch 
per minute. Three techniques were employed, each of which 


provided limited Success in critical load identification: 


ae Separation of Grid Lines 


The method most frequently used consisted of 
5X optical observation of the crack-tip grid lines. With a 
high-intensity lamp illuminating the crack-tip area, it was 
possible to determine when the crack penetrated the vertical 
line located at the crack tip aS shown in Figure 6. This 
method wcrked equally well with both liner and propellant. 
In general, the crack propagation was easier to detect in 
the propellant samples because elongation waS erratic, 
Sharp-tipped, and the material surface highly reflective, 
thereby enhancing optical techniques. The liner 
deformations at fracture levelS maintained an elliftical 


crack-tip gecmetry. 
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b. Separated Grid Limieeat eemrer edie cadr 


Figure 6. Crack Propagation Grid Lines. 
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b. Optical Observation of Crack Surface 


An alternate optical technique was used with 


the liner samples when the boundary displacement angles were 


greater than 45°, Pre-fracture deformation opened the 
initial cracks approximately 0.2 inch. The surfaces were 
painted with drawing aot when the opening reached 


approximately 0.1 inch. The high intensity light was aimed 
along the crack line at the internal crack tip. The ink 
dried quickly and the crack elongation was identified by 
parting of the blackened fracture surface, A 5X eyepiece was 
also used with this technique, but the phenomenon could 
quite well be observed by the naked eye. This technigue did 
not work witn the propellant since pre-fracture strains did 
not sufficiently open the crack to allow preparation and 


observation. 


Cz Fluctuation of Load Data 


A third technique was quite successful with 
propellant samples due to their erratic crack-elongation 
behavior. The data acquiSition system was adjusted to 
maximum gain “with zero offset so that recorded data 


emcompassed only a narrow band of force level. (Typical 


ZS 


— _ 
- © » - s 





fracture loads were in the order of 15 pounds. 
Amplification to produce 5 pounds’ full-scale aeSensler 
deflection with a 10-pound zero offset Wokemmic ec.) Tension 
up to the critical load produced a smooth data trace. 
Propagation of the crack was easily detected by fluctuations 
of the data trace. (Friction of the loading device 
typically produced load data-trace perturbations in the 
order of 0.05 pound peak-to-peak (PTP) or 0.10 inch during 
pre-~fracture loading. The propellant samples! erratic crack 
propagation produced fluctuations in the order of 0.20 pound 
Bees or 0.40 inch on the recorded trace.) This technique 


correlated well with the optical method of grid line 


separation. 
4. Trajectory Measurement 
The trajectories were optically determined 


following specimen removal from the test rig, and relaxation 
to original shape. The liner fracture was characterized by 
a well-defined straight trajectory. The propellant 
trajectories were not smooth, but zZig-zagged about a rather 
well-defined constant-direction path. Samples of each are 


shown in Figures 7 and 8. 
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IV. COMPUTER METHODS 


A. FINITE-ELEMENT SOLUTION 


1. General Description of PSELST 


The analytical solution for the elastic model of 
the labratory tests was generated by the finite element 
method. The basic program, PSELST (Plane Stress_ Elastic 
Analysis using Linear Strain Triangles) was written by Dr. 
Carlos A. Felippa [3] at the University of California 
(Berkeley, June 1966) and reprogrammed to allow initial 
displacement boundary conditions by J. P. Malone (4]j at the 


Naval Postgraduate School in July 1968. 


Displacement-compatible finite elements are used to 
ensure convergence. PSELST computes in-plane deflections 
and stresses at selected sites on the body resulting fron 
in-plane loading. The basic mesh element is a quadrilateral 
composed of four 6-nodal-point linear strain triangles, the 


center point being the centroid (Figure 9). 
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Quadrilateral . Diraene re 


(4 Corner Nodes, 8 External Nodes) (35. Corner Nodes 
(#9 = Centroid) 6 External Nodes) 


Figure 9. Linear-Strain Triangles of Basic Mesh Elements. 


A detailed description of the program may be found in 


Reference 4. 


Ze Grid Mesh Design 


— — SS == 


Generation of the finite-element grid mesh is the 
Single mest important (and time-consuming) aspect of this 
solution method. While certain ground rules are established 
for efficient grid design, selection of an optimum pattern, 


1f such exists, requires considerable trial and error. In 
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this study twelve patterns were examined which were 
essentially element-size variations of four basic designs: 
1) Rectangular elements. 
2) Rectangular elements with triangular elements 
in the ViCiINleyeO feaumne crack tip. 
3) Quadrilateral elements with polar symmetry 
and triangular elements at the crack tip. 
4) Quadrilateral elements with polar symmetry 
and finite radius at the crack tip. 


Schematics of each are shown in Figure 10. 


B. CALCULATION OF INTERNAL STRAIN ENERGY 


The internal strain energy may be computed by 
integration cf the element stresses or displacements or, in 
the case of a conservative systen, by computing the work at 
the structure boundaries. The latter technique requires 
considerably less” effort, and since poundary stresses and 
displacements were provided by PSELST, a numerical 
integration = aeaeed uSing Simpson's method was employed to 
compute the strain energy. The change in strain energy 
resulting from crack Srensarwon was computed by repeating 


the entire solution process with an extended initial crack 


under identical loading conditions. 


30 








le 


=. 





ome Polar Quadrilaterals with Finite Crack Tip Radius 


Bieure 10. Mesh Grid) Patterns ate Grace lip. 


Sd 





Ve “fNALYSioennD RESUGES 


» 


A. EVALUATION OF FINITE-ELEMENT METHOD 


Evaluation of the finite-element method as a solution 
technique consisted essentially of validating the solution 


in two basic areas: 


ii General Beiavton Of thes so Litton 


Similar to the stress solutions, the displacement 


fields of Williams are a product solution; that is, 


& 
i 


a ee (9a) 


< 
tl 


£ e 9b 
2") 9, (6) (9D) 


where f(r) and g(8) are given in Equations (2) and (4). He 
may check the : and @ behavior independently by plotting the 
computer Sommen one in the. Vieewnityeror. ~-thevcrachet:lpeat 
constant 8 or r. If agreement with the Sigel call solution 
exists, this process also establishes the bounds of validity 


for the classical solution as represented by the singular 
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term of the series. A typical set of displacement curves 
for (v) at ®@ = constant is presented in Figure 11. Similar 
results were obtained for the displacements (u). This shows 
excellent correlation for a one-inch crack within a radius 
of 0.004 to 0.100 inch from the crack tip for the gecmetry 
shown. A representation of the O0-dependence for a boundary 
displacement angle of Y = 37.5° is shown in Figure 12. This 
also correlates well with the classical solution. The data 
shown in these figures represent solutions from the grid 
pattern of Figure 10c (polar quadrilaterals with triangles 
meee crack tip). This pattern produced the best solution 
agreement with the classical in the vicinity of the crack 
tip. Determining the optimum grid pattern was essentially a 
task of evaluating the solution convergence to classical fC 


and § behavior. 


ae Simultaneous Solution of Classical Equations 


A correct soluticn requires = and ae to be 
constants for any given loading of the fixed geometry. It 
is therefcre possible to solve for these quantities by 


simultaneous solution of two linear algebraic equations 


fo 





v - Displacement (inches) 





re Distance rome Gea G ke somnelieiic..) 


Figure 11. Typical Displacement Curves at 9 = Constant. 
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Figure 12. Typical Displacement Curves at r = Constant. . 
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expressing the stresses or displacements at two points 


Within the field. In the region of classical behavior, the 


values oktained should be independent of the points 
selected. Since the solution by this method is dependent 
upon adherence to the classical formulation, accurate 


results rely upon upon precise rc and @ behavior. 


Db. PNCEG Ye Relecaseanaeces 


A second techniyue for evaluating K and K 

ie TLAE 
utilizes the energy relationships of Irwin. The work at 
the moving boundary is santnees for two solutions under 
identical loading conditions, one having an incremental 
crack extension. Since meSh grid generation 1S complex and 
time consuming, and relocation of the nodal point 
coordinates for a polar grid is undesirable, a simple schene 
was used to provide crack extension: the nodal point 
coordinates of the free ends were displaced a Small amount, 


thereby effecting a change in crack length with minimum grid 


modification. This iS Shown schematically in Figure 13. 


Because the solution is linear, and K_ and — depend 

I 
only on their respective displacement components, it is 
sufficient to evaluate the energy release rates at the two 


orthogonal cases of pure opening and pure Shear. 


35 


=, i = = 





> 7 7 , &% @ y- 





eeee 





ere | | | fol. 


BOA Se \ 
oe 


b. Extended Crack Grid Pattern 
mememiemts 1, 2, 3, 4 stretched; Elements 5 and 6 compressed) 
(Relocate 6 Nodes at left end, 3 at right end) 





Freure (3. Scheme for Extendéd-Grack Neca 
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This technique is preferred since it is, in effect, an 
weverall’ or "average”™ solution and is not Sensitive to 


local irregularities within the field . 


om Simplified Analysis of an Infinite Strip 


A third method of evaluating K and kK 

i ame 
employs the simplified analysis of the infinite strip as 
developed by Lindsey [5]. His solution assumes that a 
semi-infinite crack iS inserted in an infinite strip and the 
upper and lower boundary stresses decay to zero at some 


finite distance to the left of the crack tip as shown in 


Figure ao. 





Pagure 14. Infinite Strip with a Semi-intinite Crack- 


- 


"To the far left, the sheet is completely 
unloaded and stress free. To the far right it is 


under constant stress as if there were no crack, 
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and in the center there is a transition region of 
unkncewn Character between the two constant fieids. 
If the crack advances by an amount Ac, the 
transition region moves ahead Ac. The constant 
stress field is reduced in dimension by Ac and the 
unstressed region is increased by the same amount. 
The strain energy lost fron the constant stress 


field by the crack moving ahead Ac is 


AU = 40yey*Z2btAc 


where t is the thickness of the sheet." [5] 


In the limit of infinitesimal increments in 


length, 


~ 


(10) 


crack 


(11) 


Imposing displacement boundary conditions, Lindsey 


shows that 


vie) oe 





and in a fparallel development for K 


Il 


Uo 


E 
K = a 
= Y2(l+v) b 


(12) 


(13) 





where v. andu are the initial displacements. Although the 
O O 
geometry studied herein is far from the infinite-strif case 


of Lindsey, his results differ by less than two percent 


from the finite-element energy solution in the case of K. 
I 


For the parallel case of ee the agreement 1s not’ so 
close, the finite-element solution to oe being ten percent 
dower than that of the France Strap. The disagreement 
may be accounted for by examining the assumption of a 
decayed stress field to the left of the crack tip. Ins che 
opening mode, the normal boundary stress did in fact decay 
to zero within a short distance of the crack tif. The 
finite-element stress solution for this case is shown in 
Figure 15a. Extending the crack simply translated the 
gradient region and, we effect, the result was a 
foreshortening of the constant stress field by the amount Ac 
as assumed. 

In the Mode II case, a boundary Shear distribution of 
the PSELST solution is shown in Figure 15b. It can be seen 
that the shear stress does not decay to zero for the _ short 


geometry studied, but rather remains at a significant level 
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a. Normal Stress in Mode I Loading. 
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b. Shear Stress in Mode II Loading. 


Figure 15. Stress Distributions on Specimen Upper Boundary. 


even to the specimen's edge. Elongating the crack therefore 
has two effects: the uniform stress region is reduced by an 
amount Ac as in the case of Mode I; and, the gradient is 


extended by the amount Ac. Therefore, the resultant 
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boundary work, which equals the change of internal strain 
energy, May be considered as the reduction of the uniforn 
stress over a length Ac plus that due to the redistribution 


of the shear gradient to the left of the uniform region. 


Studies indicate, for the geometry considered herein, 
that this redistribution of shear stresses accounts for a 
ten percent reduction in the value of K from that of the 
simple case wherein the distribution translates but does 


not change. 


B. DETERNINATION OF THE LOADING ANGLE 


Except in the pure opening case, the resultant forces 
on the upper and lower specimen boundaries are not 
co-linear. This couple must be balanced by side loads. The 
finite-element solution is used to determine the magnitude 
of these unknown forces by establishing the distribution of 
the boundary stresses. From Figure 16 the relaticnship 
BS on the force resultant component along the direction of 
displacement (test data) and the orthogonal components 
corresponding to Mode I and Mode II can be determined as 


follows: 


41 






ee 


Vee oo ©. 


fy - 


-_ 


Tag 


are - 





Figure 16. Diagram of Resultant Forces along Specimen 


Boundary. 
RF, = Bsj - meas i = C1 F 14 
var sina =. Sina = Uj) Fneas (14a) 
cos (a-y) 
F Fcosa ‘meas cosa CF (14b) 
x = =_eoe}-<----:er-—— > = 
: cos(a-y) 2*meas 


If the specimen were not rigidly mounted in the 
fixture, the resultant force would parallel the 
displacement and the factors a and = would equal Sina and 
cosa, respectively. Figure 17 illustrates the distortion of 
C and a from the preferred case of sine and cosine 


functions. The relationship between boundary displacement 
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angle and the resultant force loading angle is shown in 


Figure 18. 


Force Resolution Factor 





Oe 20 40° 60° 80° 


y - Displacement Angle (Degrees) 


Figure 17. Force Resolution Factors. 
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y - Displacement Angle (Degrees) 





Oe 20° 40° 60° 80° 


a - Force Loading Angle (Degrees) 


Figure 18. Displacement Loading Angle vs. Force Loading 


Angle. 
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Cs EXPERIMENTAL EVALUATION OF STRESS INTENSITY FACTORS 


For a linear system, the strain energy is proportional 
to the displacement squared. In this case it is useful to 
consider the separate energy components due to pure opening 
Or pure shear in order to establish a Se envelope. From 


Irwin we have: 


5 aeenre 

A parallel analysis of the Shearing mode, which assunes 
the crack extends an infinitesimal distance without 
changing direction, yields the same expression for = 
Using the finite-element solution for the strain energy 
release rate, AU/AC, the following relationships are 


established: 


(16a) 


E 
Ky 7 42 2S ee 
i (T96) Yo 


(16b) 


jz 
S62 47 eee 
KrT Fae! 26 


where v.andu are the displacements in the Mode I and Mode 
oO te) 
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Ii directions respectively. For the linear system, the 
force components are proportional to the displacements. The 
proportionality constants are determined by the 
finite-element solution. For the geometry considered, we 


find that: 


x ,196 
qa een 17 
Yo = 3q50q ? ne 
Fy 196 
Yo ~ Toaq ee ? ee 
Equations (16) then become: 
Ky = -1400F), (18a) 
KT) = eco (18b) 


Using the experimental load data at fracture, the 
stress intensity factors of single material samples were 


computed and are presented in Figure 19. 
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b. Single Material, Offset Crack. 


Figure 19. Single-Material Fracture Envelopes. 
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HEdoOgin and Sih foi, studyvings the £Lracturesorea brittle 
sheet (plexiglass) with oblique central cracks, found the 
failure envelope generally elliptic in ake space. 


I If 


However, K (critical value of kK where K = 0) was 
EC af Shae 

Slightly larger than K (critical value of K_), which is 
ee If 


Opposite to the viscoelastic results of this study. 


D. PAXINUM PREINCEPAL STRESS DIRECTION 


Three basic technigues may be applied to determine the 
directicn of maximum principal stress. The simplest and 
most obvious is to solve for the angle using the computed 


stress values at the crack tip nodal point. 


Depending on the grid pattern, from four to twelve 
elements may be conmnon to the tip node. A 
crack-tip-centered polar coordinate grid with 30° segments 
was used. Increasing the number of elements resulted in 
elongated elements in the vicinity of the crack tip. For a 
well-behaved solution, included angles of less than 30° are 
not recommended in order to limit the element aspect ratios 


to a maximum of 2:1. 


On the cther hand, if the polar segment angles are 
large and the quadrilaterals are kept low in aspect ratio, 


as is desired, the element size will grow too rapidly, as fC 
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increases, to allow a sufficient number of elements in the 
Preinity Of the crack tip for solution eee The 309 
segments used herein represent a compromise which allowed 
Sufficient element quantity in the crack tip region while 
Maintaining a reasonable aspect ratio. It is reccgnized 
that the solution encompassed by the inner ring of elements 
1S questionable because of the crack tip singularity. It 
has been cbserved that even a well-conditioned element grid 
will not yield solution convergence Within the 


Singularity-adjacent elements. 


It 1s interesting to note that, in spite of the 
solution singularity at the sharf-tipped crack node, the 
principal stress direction calculated from the averaged 
stresses at that node coincided with the test data and the 
blunt-tip numerical solutions. The averaged stresses were 
computed from the solutions of all elements common to the 
crack-tip node. An example is shown in Figure 20. For a 
boundary displacement angle of 309, the classical sharp~tip 
naxinum prameci pal stress direction is 50°, The 
finite-element Sharp-tip solution approached this value as r 
decreased. However, the solution diverged within two 


elements of the crack tip to a value of 339, which agreed 
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Begure 20. Comparison of Sharp andeelunt-lf1p Solutrons. 
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With the experimental trajectory for this loading. The 
blunt-tip solution approaches the experimental result, but 
it also diverged within the last element to a value of 49°. 
In this case, the singularity is absent and the solution 
divergence is typical of finite-element solutions at a free 
surface. At the crack tip, the finite-element solution 
cannot match the Singularity of infinite stress. The values 
obtained, although fictitious, represent a solution limit in 
any radial direction aS r approaches zero. It is quite 
possible that the average of these may be considered a 
"smearing" of the solution in the vicinity of the crack tip 
and that the direction of maximum principal stress can be 


calculated from these average stresses. This result was 


observed to be independent of the grid pattern. Time 
limitations precluded further investigation of this 
eumlOosSsity. 


A second technique to determine the maximum principal 
stress direction utilizes the finite-element solution in the 
crack-tip near-field, excluding the singularity region. 
Once the limits of clasSical theory behavior have been 
established, the field within this region may be examined at 
fixed radius to determine the maximum stress locaticn and 


its corresponding direction. If, in fact, the solution is 
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Classically well-behaved, the angular location of the 
Maximum stress elements should be constant althcugh the 
stress magnitude changes inversely with the square root of 


the radius. 


Superposition of Eguations (1) and (3) for the case of 


combined loading, yields 


oy, = £, (r) +g, (6,K,,Ky7) (19a) 
oO, = £,(r)+g7(8,K,,Ky7) (19b) 

- A 19¢@ 
Oxy f(r) g,(6,K,,K,) ( ) 


where orientation of the element iS Shown in Figure 21. 


For the princifal stress direction: 


20 IE OK 0§ 
See arctan) - barctan (see) (20) 
Oy “Ty fo°82°41°81 
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Figure 21. Princifle Stress Orientation Near Crack Tif. 
Samce f£f (cr) = £ (cr) = £ (r 
1 2 jue 


Zoe ye ie 
é = sarctan{———-—_} (21) 


which 1S invariant with respect to radius. The orientation 
is unchanging along any radius, the maximum stress having 
been determined at some distance from the crack tip. In the 
limit aS xr approaches zero, even though the finite-element 
solution diverges at the singularity, one may interpret this 
orientation as being the maximum princirfal stress at the 
crack tip and hence predict the trajectory for a propagating 


Sharp-tipped crack. 
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A third method involves an analytical solution to 
Williams’ equations. In the direction of maximum principal 
stress, the shear stress is zero. From Equation (7c) we 


then have, in the direction of maximum principle stress, 


K,sin8 + Koc os hs) = 0 (22) 


where Bis the predicted trajectory angle. Solution of this 
is dependent only on the ratio of K tO K which was 
I IT 


experimentally determined by fracture testing. The angle 8 


corresponds to the sSharp-tipped crack solution. 


Se EFFECTS OF CRACK TIPsRAviUs 


¢ 


Because of the large strain levels, the geometry cf the 
Originally sharp-tipped cracks had deformed sufficiently to 
approximate an elliptic shape at fracture. Figure 22 shows a 
typical sample wherein the propagating crack-tip geometric 


deformaticn can be seen. 


The character of the crack-tip solution is dependent 
upon its relative bluntness. Since the crack tip radii of 
curvature became significant at fracture, the finite-element 
solution was adjusted accordingly to model this effect. It 


was determined that the direction of maximum principal 
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Stress (predictea trajectory direction) was relatively 
insensitive to tip radius within a range of 0.004 to 0.100 
nec hs (The minimum radius examined was 0.004 incn). Figure 
23 illustrates the sensitivity otf predicted trajectory angle 


to ciack tip radius for three selected load combinations. 





Figure 22. Crack-Tip Deformation at Fracture Lcads. 


Examination of the liner samples during fracture 
indicated tip radii in the order of 0.006 to 0.010 inch. 
Since the finite element solutions demonstrated trajectory 
insensitivity within this span of radii, the experimental 


results were compared with the blunt-tipped crack solutions. 
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Figure 23. Effect of Crack-Typ. Radius on, fra Jeerory. 


The results cf Single~material Samples, which show excellent 


correlation, are presented in Figures 24 and 25. 


Growth of the tip radius was not observed in the 
propellant samples during the fracture tests due to the 
erratic nature of the crack propagation. Unlike the liner, 
the propellant fracture tip would) "Streteh", and) then 
suddenly elongate a small distance in an interrupted 
sequence. The resultant trajectories, however, were 
identical to those of the liner under Similar loadings, and 


indicate that an effective blunting of the crack tip may 


have occurred. 
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Figure 25. Single-Material, Offset-Crack Trajectories. 
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Es BE-HATERLAL RESULIS 


Combined loading of the bi-material samples produced 
the fracture envelopes shown in Figures 26 and 27. The 


fracture trajectories are presented in Figures 28 and 29. 


The propagating cracks did not, under any load 
orientaticn, continue across the liner-propellant interface. 
Trajectories that ran to the interface abruptly changed 
direction and continued along the interface. Cracks 
Originating at the interface, or having progressed to it, 
would not depart from the interface regardless of Ilcading. 
Post-test examination of the samples revealed that the 
interface cracks propagated within the propellant; that is, 
propellant was observed on both sides of the fracture 


surface. An example is shown in Figure 30. - 


Solution discontinuities at the interface region of 
bi-material finite-element studies were not sufficiently 


resolved to allow trajectory prediction for these samples. 
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beoure 26. bBi-haterial, Centertive-Crack Fracture Eivertope. 
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Figure 28. Bi=eMacerial, Cfiscet -craeweitajectories. 
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Figure 29. Bi-Material, Centerline-Crack Trajectories. 
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Figure 30. Photograph of Bi-Material FPracture. 
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Vs CONCLUSIONS 


The application of PSELST to crack tip analyses has 
been shcwn to be a highly effective solution technique for 
Single material structures. Although the stress field 
Gamnot match the classical theory at the singularity foint, 
the soluticn converges rapidly and agreement with the 
Classical solution occurs within one or two elements fron 
the crack tif. In addition, solution of the stress field by 
this program establishes the range within which classical 


behavior may be expected. 


Because of the solution technigue, the stress 
discontinuity at bi-material interfaces limits the 
effectiveness of PSELST for stress analyses in the immediate 


~~ 


vicinity of the interface. 


It has been shown that viscoelastic trajectories may be 
predicted using linearized analyses of elastic materials to 
determine the direction of maximum principal stress. In the 
Mretnity of the crack tip, the directicn of naxinum 
principal stress is invariant with respect to radius and, in 
the case cf single materials, independent of the elastic 


modulus. 
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The large deformations associated with viscoelastic 
fracture result in an actual or effective "blunting" of the 
crack tip. Corresponding analyses must account for this 
finite radius since the stress field does not coincide with 
the sharf-tip solution. It has been found that the 
direction of maximum principal stress iS invariant with tip 
radius within a range of 0.004 to 0.010 inch for a 1.0-inch 


crack length. 


63 





7 


ee eee ee ees An 


VI. RECOMMENDATIONS FOR FURTHER STUDY 


Several areas of related interest have keen identified 
mabang this study as having potential consideration for 
further investigation. Those which may complement this and 
subsequent investigations inciude: 

A refinement of the viscoelastic sample fabrication 
process to eliminate the "soft" regions of incomplete cure. 
Uniform samples would enhance identification of the material 
properties, since strain data could easily be acquired 
Simultaneously with the force data during fracture tests. 

Sample testing should be performed in a limited time 
Span, with environmental control during storage to minimize 
the aging and absorption of moisture, which greatly affect 
the material modulii. 

The resultant critical loads should be studied to 
determine ccmpliance (if any) with the known fracture 
theories. 

The test fixture, which was rigidly mounted to the test 
machine, should be redesigned to incorporate a “TEloauamng= 
attachment which cannot transmit moments, thereby causing 
the resultant force to parallel the displacement and hence 


equal the measured force. 
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The angular range of combined loadings should be 
expanded. This study ranged from displacement loading 
angles of 0° to 90° but further studies should include the 
entire range of 360° if possible. 

The combined load tests should be evaluated with an 
elastic sample of known material properties. This would 
allow a rapid and accurate evaluation of the test procedures 
and analyses. 

The agreement between the blunt-tipped crack trajectory 
angles and the sharp-tipped solutions, determined by 
averaging the stresses at the tip from all elements common 
woelt, Warrants further study. 

An expanded analysis should be performed to investigate 
the effects of large displacements in the vicinity cf the 
Seack tif. 

The finite-element program should be modified to 
include special elements, such as long thin elements and 
crack-tip elements, to enhance fracture study. 

The refinement of technique to minimize solution 
discontinuities at interfaces of different materials should 
be studied. 

A finite-element analysis of the offset crack samples 
with a finite radius crack tip should be conducted to 


complement this study. 
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